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ABSTRACT 

Context. The A-variance analysis, introduced as a wavelet-based measure for the statistical scaling of structures in astronomical maps, 
has proven to be an efficient and accurate method of characterising the power spectrum of interstellar turbulence. It has been applied 
to observed molecular cloud maps and corresponding simulated maps generated from turbulent cloud models. The implementation 
presently in use, however, has several shortcomings. It does not take into account the different degree of uncertainty of map values 
for different points in the map, its computation by convolution in spatial coordinates is very time-consuming, and the selection of the 
wavelet is somewhat arbitrary and does not provide an exact value for the scales traced. 

Aims. We propose and test an improved A-variance algorithm for two-dimensional data sets, which is applicable to maps with variable 
error bars and which can be quickly computed in Fourier space. We calibrate the spatial resolution of the A-variance spectra. 
Methods. The new A-variance algorithm is based on an appropriate filtering of the data in Fourier space. It uses a supplementary 
significance function by which each data point is weighted. This allows us to distinguish the inffuence of variable noise from the actual 
small-scale structure in the maps and it helps for dealing with the boundary problem in non-periodic and/or irregularly bounded maps. 
Applying the method to artificial maps with variable noise shows that we can extend the dynamic range for a reliable determination 
of the spectral index considerably. We try several wavelets and test their spatial sensitivity using artificial maps with well known 
structure sizes. Performing the convolution in Fourier space provides a major speed-up of the analysis. 

Results. It turns out that different wavelets show different strengths with respect to detecting characteristic structures and spectral 
indices, i.e. different aspects of map structures. As a reasonable universal compromise for the optimum A-variance filter, we propose 
the Mexican-hat filter with a ratio between the diameters of the core and the annulus of 1.5. When the main focus lies on measuring 
the spectral index, the French-hat filter with a diameter ratio of about 2.3 is also suitable. In paper II we exploit the strength of the 
new method by applying it to different astronomical data. 

Key words. Methods: data analysis ~ Methods: statistical - ISM: clouds - ISM: structure 



^ 1. Introduction promising first approach to a parametric description is given by 

exponents of scaling relations. 

The interstellar medium is highly turbulent and turbulent mo- Motivated by the similarity of observed interstellar cloud im- 

tions determine the evolution of interstellar clouds. The turbulent ages with th e structure of fractional Brownian motion (fBm, see 

pressure is partially ab le to support them against gravitational Sect. l2XTl i fractals, which are characterise d by the single num- 

collapse (iKlessen et alj |2000); turbulent shocks create and dis- ber of the exponent of the power spectrum, IStutzki et alj ([T998I) 

solve dense clumps in mo l ecular clouds or even whole clouds developed the A-variance analysis as a tool to measure the struc- 

dBallesteros-Paredes et alt |1999|), the turbulent mass trans- tural scaling behaviour of observed images. 

port modifies their chemical evolution ("Decamp & Le Bourlotl n,, » . . j- , , j- ■ 

| !,nn%K J *u • 1 * u 1 » » » A ^ ' ■ »u ■ The A-variance is a type or averaged wavelet transform that 

120021) . and the iu'egular turbulent structure determines their pen- , . .-'^ ^ ■ , , 

. u J- irr- 1- — i T~^\i'^^\^\n^ t-u mcasurcs the Variance in a structure /(r) on a given scale / 

etration bv UV radiation (IZielinskv et al.L I2000). Thus, the com- , r-, ■ ■ , ■ • ■ , , ^ 

, J • , , — 71 1 ^ 1,- £ ^ 1. 1 by filtering it by a spherically symmetric down-up-down func- 

plex dynamic structure on all scales resulting from turbulence .■' ^ ■ , ,nr- , or, /- ^r^r^r^-. t,, . . , 

I. ■ ^ ^ ■ T ^- c . f.u . u • tion of size / (Zielinsky & Stutzki 1999). The A-variance analy- 

has important implications for many aspects of the astrophysics . ,-„ ,., , , . , , 

c^u SIS was successfully applied to several observational data sets: 

of the interstellar matter. „ -rr n i „ i , .c i ^ ^ i 

Stu tzki et all n99m studied a CO map of the Outer Galaxy, 

Whereas many observations reveal the complexity of the iBensch et alj (1200 lb investigated a series of nearby star-forming 

structure of the interstellar medium, most models of interstel- clouds and a number o f nested maps in different CO isotopes 

lar clouds are still based on simple geometrical configurations, from the Polaris Flare, iHubeii (12002) perfor med a systematic 

A first step towards a better understanding of interstellar turbu- study of a large set of Galactic CO maps, and Sun et al.l (l2006h 

lence and towards building more realistic models of interstellar analysed maps of the Perseus molecular cloud taken in var- 

clouds is to identify model structures characterised by a limited ious tracers and including the analysis of velocity channels. 

set of parameters which can be quantified by comparison with The intensity maps of most clouds resulted in power-law A- 

observed cloud images. As many aspects of o bserved in terstellar variance spectra with exponents between 0.5 and 1.3. Mac 

clouds can be described bv fractal properties (ICombesll2000i) . a Low & Ossenkopf (2000) and Ossenkopf (2002) applied the 
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A-variance analysis to simulations of interstellar turbulence to 
compare the scaling behaviour of the simulations with that of 
observed maps. It became however obvious that, aside from the 
spectral index, deviations from a power law on particular scales 
should be studied as well because they provide significant in- 
formation on the physical processes on these scales. Thus the 
A-variance analysis is to be optimised with respect to its capa- 
bilities of the corresponding scale detection. 

We propose in this paper a number of improvements to the 
A-variance optimising its sensitivity, its applicability to arbitrary 
data sets, and the speed of its computation. The critical quan- 
tity for the detection of pronounced scales in a structure is the 
shape of the wavelet filter function. T he spherically symm etric 
down-up-down function introduced bv lStutzki et alJd 19981) is an 
obvious first choice. However, other wavelet shapes offer attrac- 
tive alternatives. 

For infinitely extended or for periodic structures the fastest 
way of numerically calculating the A-variance is given by a 
Fourier transform of the image. However, observed maps typi- 
cally have a finite size, often even cutting the observed clouds 
at the map boundary, and Fourier-based methods run into the 
well known problems of artificial s t ructu re being introduced 
by these edge effects. iBensch et alj (12001,) thus implemented 
the A-variance by a numerical treatment in the spatial domain. 
Calculating a two-dimensional convolution in the spatial do- 
main, however, results in a rather slow computation. An addi- 
tional complication in observed data comes from the fact that 
the signal-to-noise ratio is often not uniform across the mapped 
area. A particular example of maps with strongly variable data 
reliability are line centroid velocity maps. Here, the accuracy 
of the centroid velocity alway s depends on the line intensity. 
iMac Low & Ossenkopf! ( l2000l) have shown that the "traditional 
A-variance analysis may fail in this case. To relieve these prob- 
lems and concerns, we introduce a supplementary function into 
the A-variance analysis which is used to weight the data points 
in the spatial map according to their significance. This helps to 
derive correct contributions of data points with a different signal- 
to-noise ratio to the structure information on a particular spatial 
scale and it allows us to calculate the A-variance in Fourier space 
and thus to make use of the numerical advantages of the fast 
Fourier transform algorithm. 

After revising the fundamental properties of the A-variance 
and defining appropriate images to test the method in Sect. 2, 
we introduce the concepts of the improved A-variance includ- 
ing a weighting function in Sect. 3 and optimise it with respect 
to the wavelet filter function in Sect. 4, where we also verify 
its performance by extensively testing it against the test struc- 
tures. A combined test of the optimised filter function and the 
significance function in the case of noisy data is presented in 
the Appendix. Sect. 5 summarises our findings providing rec- 
ommendations for the optimum method and wavelet to use. In 
a second paper, we test the capability of the new method ap- 
plying it to simulations of interstellar turbulence and observed 
molecular line maps exploiting the improved sensitivity to de- 
rive general properties of interstellar turbulence. 



2. The starting point 

2.1. The A-variance 

The A-v ariance analysis was comprehensively introduced by 
ISfiitzki et al. ( 1998) and Bensch et al. (2001). Here, we only re- 
peat those equations which are essential to understand the exten- 
sions proposed in Sects. [3]and|4l Although the A-variance can 



be used in principle for an arbitrary number of dimensions we 
restrict ourselves to the two-dimensional case, i.e. the analysis 
of maps or images. 

The A-variance measures the amount of structure on a given 
scale I in a map /(r) by filtering the map with a spherically sym- 
metric down-up-down function of size / (French-hat filter) and 
computing the variance of the thus filtered map. It is given by 



(1) 



where, the average is taken over the area of the map, the sym- 
bol * stands for a convolution, and Q; describes the French-hat 
function defined as 




|r| < 1/2 

1/2 < \r\ < 31/2 

\r\ > 31/2 



(2) 



Thus the filter consists of a positive core and a negative annulus 
where the width of the annulus agrees with the diameter of the 
core and the absolute values in each of them are inversely pro- 
portional to their areas so that they both have an integral weight 
of unity Q. 

In a more general picture one can consider the filter function 
as a wavelet composed of a negative and a positive part both 
normalised to integral values of unity so that the overall filter has 
a vanishing integral. Using an arbitrary diameter ratio between 
the annulus and the core v we can write 



(3) 



/,coi"e 



/,ann 



l/(v2- 1) : l/2<\r\<vxl/2 



: |r| < 1/2, |r| > V X 1/2 



(4) 



The "traditional" A-variance filter is reproduced for a diameter 
ratio V = 3. We come back to discussing the choice of v in Sect. 
4. 

Because the average distance between two points in the core 
and the annulus of the filter is close to the length / (see Sect. 14. 2b . 
the convolved map only retains variations on that scale whereas 
variations on smaller and larger scales are suppressed. The A- 
variance as the variance of the convolved map thus measures 
the amount of structural variation on the scale I. Plotting the A- 
variance as a function of the filter size I then provides a spectrum 
showing the relative amount of structure in a given map as a 
function of the structure size. 

The filter convolution and computation of the A-variance can 
be easily performed in Fourier space where they are reduced to a 
simple multiplication and integration. This directly relates the A- 
variance to the power spectrum. If P(|^|) is the radially averaged 
power spectrum of the structure /(r), the A-variance is given by 



2n 



p(i^i) Q^im 



\k\ d\k\ 



(5) 



where is the Fourier transform of the filter function with the 
size I and k denotes the spatial frequency or wavenumber If the 



' Following the original definition by IStutzki et alJ l ll998h our A- 
va riance is higher by t he constant factor of 2;r than the definition used 
bv lBensche"taLld2001l) . 
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power spectrum is given by a power law, P(|^|) oc |^| the A- 
variance also follows a power la w cr^ oc l" with a = ^ - 2 within 
the exponential range < ^ < 6 dStutzki et aUll998lPI . 

Thus, the A-variance shows in principle only information 
that is also contained in the power spectrum. The main advan- 
tage of the A-variance method compared to the direct computa- 
tion of the power spectrum results from the smooth filter shape 
which provides a very robust way for an angular average, insen- 
sitivity to singular variations, and independence of gridding and 
finite map size effects. It provides a good separation of different 
effects based on their characteristic scale, e.g. a clear distinction 
between observational noise, structure blurring by the finite tele- 
scope beam, and the internal scaling of the astrophysical source. 
A detailed computation of the infl uence of finite map s izes and 
telescope blurring was provided bv lBensch et alj (1200 ll) . 

However, when applying Eq. (|5]l to compute the A-variance 
it inherits a main drawback from the power spectrum - it im- 
plicitly assumes a periodic continuation in the Fourier transform 
although most astrophysical observations show no periodicity. 
Sect. [3] thus deals with the implications of this assumption and 
possible ways to overcome the resulting limitations. 

2.2. Test data sets 

In order to test how the A-variance reproduces specific struc- 
tural characteristics, we have constructed a series of artificial 
data sets with known characteristics. They were either chosen 
to reproduce the typical self-similar scaling behaviour measured 
in many astrophysical observations (Combes, 2000) or to con- 
tain pronounced artificial structures with a well known size scale 
which should be clearly detected in the A-variance spectrum. 
All test data sets were generated with an intrinsic resolution of 
128 X 128 pixels. 

2.2.1. Periodic data 

As a simple test bed where boundary effects play no role we 
started with three types of periodic maps. The first type is pro- 
vided by the frac tal structures of frac tional Brownian motion 
(fBm) as used by iBensch et alj ( 1200 lb . fBm structures are de- 
fined by a power-law power spectrum P(|^|) oc and random 
phases in Fourier space. We created fBms by the Fourier trans- 
form of a Hermitian field with amplitudes following the power 
spectral index f and random phases. This procedure guarantees 
real values and periodic maps. The structure analysis in terms 
of the A-variance should recover the power spectral index of the 
fBm structure measuring a slope a - ( - 2 for this data set. 

As periodic structure with a pronounced size scale we use 
chess board like patterns where the number of fields on the board 
is varied to change the size of the structures in the map relative to 
the map size. Because the single chess fields are the only struc- 
ture in this data set, their size should appear as prominent peak 
in the A-variance spectrum. This data set gives a sharp defini- 
tion of the characteristic length scale in the spatial domain but 
contains a high contribution of high frequency modes in the spa- 
tial frequency domain due to the sharp edges of each field. Thus 
we have added as a third type of test maps data fields provided 
by a single Fourier component in both directions, i.e. the super- 
position of two orthogonal sine waves with the same periocQ, 

^ Please, note the difference to the often used energy spectrum P{k) 
which is obtained by angular integration of P{k) and has the spectral 
index = f - 1. 

^ All maps are normalised to have a side length of 1 here. 




0.2- - -0.2 

0.0 1 , , , \ , , , \ , , , \ , , , \ , , , \ \ 

0,0 0.2 0-4 0.6 0.8 

Fig. 2. Example of a filled circle structure used to test the A- 
variance algorithm. The diameter of the circle is 0.25. 



/(r) - cos(27rA:x) + sin(27r^y). In the case of a wavenumber 
k- Hi can be regarded as an fBm with ^ — > oo. The scale of the 
characteristic variation should be clearly detected but in contrast 
to the power spectrum we do not expect a single sharp maximum 
in the A-variance spectrum, because the Fourier transform of the 
A-variance filter is a Bessel function with pronounced side lobes. 

The test data sets are all characterised by one free parame- 
ter. For the fBms this is the spectral index (3. For the chess board 
pattern and the sine wave field it is the size of the characteristic 
structure or the dominant wavenumber, respectively. Examples 
of the test data sets are shown in Fig.[T] The fBm structure used 
here is characterised by a spectral index f = 3.0 and a total vari- 
ance of 1 . The chess board example shows characteristic struc- 
ture lengths between and 0.24. Its average size, integrated over 
all possible angles is 0.13. The sine wave field shown here is 
characterised by a wavenumber ^ = 8 leading to a length of 0.09 
for the maximum variation. 

2.2.2. Non-periodic data 

As real astrophysical data are hardly periodic, tests for the treat- 
ment of boundary eff'ects have to be performed on non-periodic 
data. We use two types of data sets to study these effects. 

First we select subsets of l arger (periodic) fB m structures in 
the same way as introduced bv lBensch et al.l (120011) . The subsets 
cover one quarter in length, i.e. 1 / 16 in area, of the periodic fBm 
field so that they should hardly retain any information about the 
large-scale periodicity. As the subsets are randomly chosen they 
typically show sharp discontinuities at the edges. This seems to 
exclude Fourier based methods for the analysis. In Sect. |3.1| we 
show, however, that the A-variance analysis can be extended to 
account for the discontinuities. Although it is not guaranteed that 
a subset has the same spectral index as the whole fBm struc- 
ture we will judge the value of the A-variance analysis based on 
the agreement of the determined spectral index with the origi- 
nal fBm index because this approach reflects the typical obser- 
vational strategy that high resolution observations are restricted 
to small parts of a molecular cloud but they are used to derive 
the general scaling behaviour of the cloud (see e.g. the IRAM 
keyproject "Small-sca le structure of pre-star-forming regions", 
iFalgarone et 111119981) . 

As second non-periodic structure we use a filled circle on an 
otherwise empty map. The dominant scale is the size of the cir- 
cle but as the A-variance also measures the size of the "empty" 
surroundings of the circle we expect considerable contributions 
from this area to the spectrum at large lags. By adjusting the po- 
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0.0 0.2 0.4 0.5 0.8 0,0 0,2 0,4 0,6 0,8 



Fig. 1. Examples of periodic data sets used to test the A-variance algorithm. The fBm structure on the left is characterised by a 
spectral index ^ = 3.0. The chess board has a characteristic structure size of about 0.13. The sine wave field on the right consists of 
a single Fourier component at ^ = 8, the corresponding characteristic distance is 0.09. 



sition of the circle with respect to the map boundaries we can 
test the robustness of the edge treatment in the non-periodic al- 
gorithm (in a periodic treatment the A-variance is independent 
of the position of the circle). The main free parameter of this 
structure is the diameter d of the circle. The average circle scale 
is given by njAd. We have not scanned the distance to the map 
boundary as an additional free parameter, but only performed a 
few tests for circles shifted by a integer multiples of the full di- 
ameter relative to the centre. An example with a circle diameter 
of 0.25 is shown in Fig.|2l 



3. Weighting in Fourier space 

3.1. Edge treatment 

In the introduction of the A-variance by'Stutzki.etaD(11998') we 
made extensive use of the Fourier transform, applying Eq. (|5]l 
to convolve the astrophysical image with the A-variance filter 
With the Fourier transform implicitly assuming data periodicity, 
we introduce, however, steps at the edges of a map if the intensity 
does not show the same value on both sides. Because step func- 
tions contain contributions at all spatial frequenci es they distort 
the A- variance spectrum of the original structure. iBensch et alj 
(1200 lb have shown that this effect can lead to considerable errors 
for most astrophysical structures, where a periodic continuation 
is not possible. 

The solution proposed there with the POINT and PIX al- 
gorithms was a modification of the filter function when applied 
close to the boundaries of the maps. The filter is truncated so 
that it never stretches beyond the map edges. To guarantee that 
the core and the annulus are still normalised to unity for the trun- 
cated filter, both parts are multiplied with correction factors de- 
pending on the remaining filter size. This edge treatment has, 
however, the disadvantage that different parts of a map are con- 
volved with a different filter function so that the computation 
of the A-va riance by Fou r ier tra nsform via Eq. Q is no longer 
possible. In ^ Bensch et alj (1200 ll) we thus concluded that the A- 
variance should be rather computed in ordinary space. This ap- 
proach is, however, slow compared to the computation in Fourier 
space (several hours instead of a few seconds for maps contain- 
ing more than about 100x100 pixels). 

Here, we introduce a method that combines the improved 
edge treatment with a computation in Fourier space. It is fast 
and does not introduce artificial high-frequency contributions by 
periodically wrapping discontinuities at the map edges. To ob- 



tain this behaviour the method is set up to fulfil three conditions: 
i) In the convolution a fixed filter function is to be used so that 
the convolution can be computed by multiplication in Fourier 
space, ii) To avoid effects of periodic edge-wrapping, the filter 
contributions beyond the edges of the map have to be truncated, 
iii) The normalisation of the filter discussed in Sect. 12.11 has to 
be fulfilled at each point. 

These apparently conflicting requirements can be met using 
two simple ideas. The truncation of the filter is substituted by 
a zero-padding of data and the corresponding error in the filter 
normalisation is corrected by weighting factors for the two filter 
contributions computed as a function of the coordinates in the 
map. 

Instead of truncating the filter when it extends beyond the 
map edges we increase the map by zero-padding beyond the 
edges up to the maximum filter size used. Then the extended 
map can be convolved with a fixed filter function, only providing 
zero contributions from outside the original map. In this way no 
points from other periodically wrapped parts of the image may 
then fall into the filter centred at any map position. The error in 
the normalisation of the filter introduced by this substitution can 
be computed from the convolution of the filter with an auxiliary 
map w(r), which has a value 1 inside the range of valid data and 
in the zero-padded region. Because the A-variance filter (3/ 
has to fulfil the two normalisation conditions for the core and 
the annulus (see Eq.|4]i 

ZOw = ZO('-) = i (6) 

map /.core map ;,ann 

it has to be split into the positive and negative filter parts for the 
computation of the normalisation errors. The sums extend over 
the map of valid data. In total four convolution integrals need to 
be evaluated to compute the A-variance 

Gicoreir) = /padded('-) * Qir') 

/,core 

G/,ann('-) = /padded('') * 

/,ann 

W,,core(r) = w(r) * Q(r') 

/,core 

W,,a„„(r) = w(r) * Q(r') (7) 

/,ann 

where /padded stands for the map with the additional zero-padded 
boundary region. As we use fixed filter functions and all involved 
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functions vanish at the map edges, the convolution integrals can 
be easily computed in Fourier space involving a fast Fourier 
transform and a map multiplication. Although the map treated 
in this way is larger than the original map by up to a factor three 
in each direction, the Fourier transform using an algorithm like 
FFTW, which can work on arbitrary map sizes, is still consider- 
ably faster than any convolution in ordinary space. 

The full map convolved with the A-variance filter truncated 
at the map edges is then 



Fi(r) 



W,,core(r) " Wl,,„nir) 



(8) 



and 



It is only defined where the normalisation parameters W/_c 
ann ^6 both different from zero. 

In the computation of the A-variance spectrum one has to 
take into account the reduced significance of the data values in 
the convolved maps produced by the fact that the applied filter 
becomes more and more distorted relative to the optimum filter 
when it is truncated. Using the normalisation parameters of the 
truncated filters as a measure for their significance one can add a 
significance weighting to the A-variance analysis. We define the 
A-variance no longer as the variance of the convolved map but 
weight the map points by the significance of the filter applied to 
compute the value at each point when computing the variance 



with 



2map Wl^totir) 



,annV' / 



(9) 



(10) 



The sum covers the whole (extended) convolved data field. The 
definition of the significance function as the product of both nor- 
malisation factors is somewhat arbitrary but reproduces the de- 
sired behaviour that changes in the positive and negative part 
of the filter contribute equally. We have tested different powers 
of the product but found the best agreement with the theoretical 
behaviour in fBm data sets for an exponent of just unity. 

Figure [^demonstrates the effect of the edge treatment in the 
example of an fBm structure which is periodic and for a non- 
periodic sub-map from a larger fBm structure. We have used 
three different ways to compute the A-variance. First we as- 
sume that the maps are periodic, neglecting all wrap-around ef- 
fects. For the periodic fBm this is, of course, the best assump- 
tion which should reproduce exactly the properties of the power 
spectrum used to generate the fBm. It is, however, rarely useful 
when dealing with observed data as they are in general not peri- 
odic. The second approach creates periodicity by mirroring the 
map along both axes as discussed by Stutzki et al. (1998) so that 
a larger periodic map is produced and wrap-around effects can 
be neglected. This approach, however, still results in discontinu- 
ities in the first derivatives at the mirror axes. The third approach 
uses the truncated filter as described above. 

For the periodic fBm structure we find a very close agree- 
ment of the A-variance using the truncated filter with the theo- 
retical value given by a slope a - I. The mirror continuation 
results in an apparent reduction of the amount of large-scale 
structures in the map as they are partially assigned to larger 
modes only present in the map extended by mirroring. Thus the 
A-variance slope is systematically underestimated. For the non- 
periodic structure, the assumption of periodicity results in strong 
deviations from the power-law behaviour on small scales due to 
artificial high-frequency contributions from the edges. The mir- 
roring shows again an underestimate of the power spectral index 
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periodicity used, a = 0.986 
X mirrored, o; = 0.879 
A truncoted filter, a = 0.972 




0.01 



0.1 

lag 



0.01 



^< 10" 



<> periodicity assumed, a = 0.f 

X mirrored, o; = 0.894 

A truncdted filter, a=1.01 




0.1 
log 



Fig. 3. A-variance spectra computed with three different types of 
the edge treatment for one periodic fBm structure (upper plot) 
and one non-periodic subset from a fBm structure (lower plot). 
The fBms are characterised by a spectral index f = 3.0 so that 
the A-variance should show a slope a = 1.0. The measured 
slopes are indicated for each graph in the plot. 



on large scales, and only the use of the truncated filters results in 
a good reproduction of the expected value for the slope. 

In Fig. [3] we plot the results obtained just for one realization 
of an fBm and for one sub-map from a large fBm. But the ex- 
act shape of the computed A-variance spectra can depend on the 
distribution of the random phases within the fBm and it will cer- 
tainly depend on the selection of the sub-map within an fBm. 
Thus we have repeated the computation for a se t of maps in 
Fig.[4| Analogously to the statistical treatment by 'Bensc h et al.l 
(12001 1) we vary the spectral index of the fBm structures and 
chose randomly 30 different fBms or fBm sub-maps and de- 
termine their A-variance spectra. We do not display all spectra, 
but only the resulting distribution of slopes a. The Fig. shows 
the average A-variance slope and the spread of measured slopes 
as a function of the spectral index f using the three different 
edge treatments. All slopes are computed from a power fit cov- 
ering the full data range plotted in Fig. [3] An optimum treatment 
should reproduce the relation a - ( - 2 indicated by the dotted 
line in Fig.H] 

For the periodic fBms we find that the A-variance spectra 
from the truncated filter show about the same spectral indices as 
the periodic treatment. The error bars in the periodic treatment 
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Fig. 4. Average and standard deviation of the spectral indices 
of the A-variance spectra determined for a set of fBm struc- 
tures (upper plot) and for a set of sub-maps from fBm struc- 
tures (lower plot) depending on the spectral index of the fBms. 
All spectra are computed at the same power spectral indices 
( = 2.0, 2.4, 2.8, 3.2, 3.6, and 4.0, but for a better visibility, the 
standard deviations for the truncated filter and for the mirror 
continuation are displaced in the plot relative to these values. 
In the case of the periodic analysis of the fBm maps the error 
bars giving the standard deviation of the statistical samples are 
practically zero. The dotted line indicates the theoretical index 
for infinitely large fBm structures or = ^ - 2. 



are zero because the A-variance spectrum is independent of the 
exact phase distribution and thus identical for each map of the 
sample. The standard deviation of the slopes in the truncated- 
filter treatment is always about 0.08 here. The A-variance spec- 
trum underestimates the power spectral index by 0.03 at ^ = 3.6 
and by 0.08 at ^ = 4 because the theoretical value is only reached 
for infinitely large maps and the deviations grow when approach- 
ing the asymptotic limit of o- = 4 (see Stutzki et al. 1998). The A- 
variance spectrum computed for the mirror-continuation of the 
map always underestimates the spectral index by about 0.1. 

In the case of non-periodic fBm sub-structures the periodic- 
ity assumption clearly fails. With this approach the measured 
slope saturates at about a - 1.1. The use of the simple A- 
variance without filter truncation or mirror continuation provides 
wrong results at power spectral indices above about ^ = 2.9. This 




10-" 

0.0 



A truncoted filter, 1=0.108 

X mirrored, 1=0.1 12 

<> periodicity ossumed, 1=0.113 

0.1 
log 



Fig. 5. A-variance spectra for the filled-circle map where the cir- 
cle has a diameter of 1/8 of the map size and is shifted by three 
diameters from the centre of the map. The values / indicate the 
lag of the measured peak in the spectrum. 



is due to unavoidable discontinuities at the submap edges which 
result in high-frequency contributions in the periodic treatment 
resulting in too shallow A-variance spectra. In contrast, both the 
mirror-continuation and the filter truncation provide a reasonable 
measure for the actual map structure for all spectral indices. Both 
methods are hardly affected by the discontinuities at the sub-map 
edges. The mirror-continuation always underestimates the spec- 
tral index by about 0.1 (except at ^ = 4.0). The filter truncation 
method reveals the correct spectral index for ^ < 3.0 and overes- 
timates it by 0.05 at f = 4.0. Regarding the typical error bars of 
0.15, the systematic errors are, however, lower than the scatter 
between different fBm realizations with the same spectral index. 

The same tests were repeated for map sizes ranging fr om 32^ 
to 256 ^ pixels. In agreement with the studies by Benscl Tet al.l 
(1200 Ih we found no systematic changes in the A-variance slopes 
exceeding 0.03 when changing the map size but an increase of 
the error bars from 0.15 at map sizes of 128^ to 0.25 at a size 
of 32^. This higher uncertainty for smaller maps prevents any 
significant conclusion from maps spanning less than x 30 pixels. 

To demonstrate the influence of the edge treatment on the 
A-variance spectrum of an object with a pronounced size scale 
we show in Fig. |5]the spectra computed for the map contain- 
ing the filled circle with a diameter of 1/8 of the map size. The 
peak of the A-variance spectra at 0. 1 1 falls below the circle di- 
ameter of 0.125 but slightly above the average distance between 
two points on the rim of the circle which is about 0.10. The 
somewhat higher value is probably due to the contributions from 
the "empty" environment of the circle at large lags, also result- 
ing in a decay of the A-variance spectra at large lags which is 
shallower than the T^-characteristics of uncorrected structures.. 
These contributions are dispersed over a relatively wide range of 
scales corresponding to the different distances to the map bound- 
ary in a non-periodic treatment and to the distances to the next 
circle in a periodic interpretation of this structure. As these vari- 
ations are mainly assigned to lags exceeding the map size in 
the periodic treatment, the two curves for the periodic treatment 
show somewhat lower A-variance values within the map than the 
filter truncation method where the "empty" region is constrained 
by the map size. Nevertheless, the total differences between A- 
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variance spectra using the different edge treatment methods are 
relatively small, so that either method seems to be justified for 
this case. 

3.2. Generalisation for data witli varying reliability 

The concept of weighting the A-variance computation by dif- 
ferent filter significance values can be generalised to deal with 
data, where the data points in a map are as well characterised by 
a variable data reliability. This applies e.g. to maps where not all 
points are observed with the same integration time so that they 
show a different noise level. The inverse noise RMS is an indi- 
cator for the significance of the data at different points. Many 
other observational effects may lead to a similar variation in the 
data reliability across the map. As long as the reliability can be 
expressed as a significance number Wdata('') between and 1 all 
such maps may be analysed within the concept outlined here. 
The same equations as discussed in the filter truncation are to 
be applied, but the auxiliary weight map w(r) does no longer 
consist of the values 1 inside and outside of the original map. 
It rather contains the significance values Wdata('') ranging con- 
tinuously from to 1 . The weighting factors in the A-variance 
computation W^otir) then contain the integrated significance of 
the filter-convolved data at each point. 

With this generalised concept, the A-variance analysis can 
be applied to arbitrary two-dimensional data sets. They must be 
projected onto some regular grid but they do not need to con- 
tain regular boundaries as the corresponding "empty" grid points 
only have to be marked with a zero significance. Varying noise or 
other changes in the data reliability can be expressed in the sig- 
nificance function Wdata('') which has to be constructed for each 
data set. The only remaining requirement for the applicability 
of the A-variance is the sufficiently large spatial dynamic range 
in the data. The criterion of at least 30 pixels in each direction 
for reasonable error bars of the A-variance spectrum discussed 
above has to be extended in the case of a low data significance. 
In the appendix we present measurements of the dynamic range 
over which the slope of fBms can be reliably determined in the 
case of noisy data. We find as a rule of thumb that the minimum 
map size has to be increased by one over the average data sig- 
nificance. In paper II we will apply the A-variance analysis to 
observed data with irregular boundaries and a spatially varying 
significance. 



.8 r ' ' ' ' ' " 

French hat 

g _ X Mexican hat 

.4 - 




0.4 [ , , , , \ , , , , \ , , , , I , , , , 

1.0 1.5 2.0 2.5 3.0 

diameter ratio 

Fig. 6. Factor translating the filter core diameter I into the aver- 
age distance measured by the filter as a function of the diameter 
ratio between annulus and core of the filter. The average distance 
is computed by a double integral over the core and the annulus 
of the filter 



where I is the size of the filter and v is the diameter ratio between 
the annulus and the core of the filter as defined in Sect. 12.11 The 
choice of the Gaussians guarantees that the filter gives the best 
simultaneous confinement in ordinary space and in Fourier space 
not showing any side lobes in either of them. 

As the opposite extreme we have also tested a filter measur- 
ing the difference between one point in the map and all points 
displaced by the sharp distance / relative to this pixel. The corre- 
sponding A-variance then measures basically the st r ucture func- 
tion o f the map (see e.g. iMiesch & Balm Il994t iFrick et all 
2001'). However, becau se those test s only confirmed the results 
from Ossenkopf & Mac Lowl ( |2002|) and Frick et al. (2001) that 
the structure function is relatively insensitive to distortions of the 
power spectrum on particular scales we excluded this filter from 
the following studies. 

Thus, we restrict ourselves here to the Mexican- and the 
French-hat filter, where we vary in both cases the diameter ra- 
tio V. In this way we test both the influence of the general filter 
shape and of the ratio between core and annulus for each filter 
on the computed A-variance spectra. 



4. Filter optimisation 

4.1. The A-variance filter function 

All examples given above were computed with the fixed filter 
function of a French hat with a diameter ratio between the annu- 
lus and the core v = 3. An obvious question is whether we can 
improve on the A-variance by using a different diameter ratio 
and/or a different filter function. Due to its discontinuity in the 
normal space the French hat has high frequency lobes in Fourier 
space. Alternative approaches should use smoother functions in 
ordinary space to obtain a better confinement in Fourier space. 
As a smooth example we implemented a "Mexican hat" consist- 
ing of two Gaussian functions: 



/,core ^ 
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i2(v2 - 1) 
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4.2. The effective filter length 

Taking into account the finite width of the core and the annulus 
of both filter functions it is not obvious on which scale variations 
are actually measured when using a filter with size I. Here, we 
compute this scale based on the geometrical properties of the 
filter. 

In the French-hat filter we measure the average distance be- 
tween a point in the core and a point in the annulus, providing the 
average scale on which a structure variation in the map should 
be measured. For the Mexican-hat filter, the computation of the 
average distance includes the additional weighting of each dis- 
tance by the product of the positive and negative filter values. 
This reflects the effect of the convolution of a map with this fil- 
ter 

Figure |6] shows the resulting effective filter length relative 
to the filter size / as a function of the diameter ratio v for the 
French and the Mexican hat. The length scale traced by the filters 
is approximately a linear function of the diameter ratio between 
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Fig. 7. A-variance spectra computed for the sine wave field with 
a wavenumber k - ^. The fiher truncation is apphed in both 
cases for the edge treatment. The peak scale is given as I, the 
logarithm of the ratio between upper and lower lag where the 
A-variance has dropped to 1 /2 of the peak value is given as w. 

the core and the annulus of the filter. Using a least-square fit we 
obtain the coefficients 

kff _ r 0.29 V -H 0.26 for the French hat 
T ~ \0.41v-h0.46 for the Mexican hat. 

The original A-variance definition using a French hat with 
V = 3 gives an effective length of 1.12 times the core diameter 
/. Thus all scales computed previously with that filter should be 
shifted by the factor 1.12. This is only a small correction, not 
changing the conclusions in any of the papers that have used 
the A-variance so far In all plots shown in this paper, we use 
the effective length as the lag of the A-variance to allow a direct 
comparison of the spectra independent of the filter used. 

4.3. Filter evaluation 

The optimum filter to be used in the A-variance analysis has to 
fulfil two criteria: the correct detection of pronounced size scales 
in the maps and the exact determination of the scaling exponents 
of the contained structures. 

4.3.1. Scale detection 

In the detection of pronounced scales the maximum of the A- 
variance spectrum should fall onto the correct lag corresponding 
to the structure size. Moreover, the signature of the pronounced 
scale in the A-variance spectrum should be as sharp as possible 
with a high contrast relative to other scales. As test images we 
used the chess board field, the sine wave field, and the filled 
circle field. 

To illustrate the behaviour we plot the results for the sine 
wave field. Fig. [T] shows the A-variance spectra measured for 
a field with ^ = 8 by the French and the Mexican-hat filters. 
In both cases the filter truncation method is used. The dominant 
scale is detected as a peak in the A-variance spectra, the peak po- 
sition at about 0.07 is somewhat lower than the scale of the max- 
imum variation l/(ky/2) - 0.088. We found this small shift by 
about 20 % in all spectra. When using the effective filter length, 
the peak position is very constant independent of the filter type 
and its diameter ratio. 



On scales above the peak, the spectra obtained for the two fil- 
ters deviate considerably. The French-hat filter produces ripples 
at large lags. From the positions of these ripples we find that 
they reflect the side lobes of the Bessel function, representing 
the Fourier transform of the French-hat filter This must not be 
misinterpreted as a detection of large-scale structure in the map. 
The side lobes detect the single Fourier amplitude from k - ^ 
also at other effective filter sizes. This can be seen when chang- 
ing the diameter ratio v in the filter. The peak remains at the same 
position but the ripples in the A-variance spectrum move corre- 
sponding to the changes in the side lobes of the Bessel function. 
Decreasing v relative to the value of 3.0 results in a steeper de- 
cay above the peak and an increased number of ripples at large 
lags. When reducing the diameter ratio v for the Mexican-hat fil- 
ter the decay above the peak also steepens and we obtain some 
flattening at the largest lags in the map. In general we can either 
achieve the sharp peak and the artificial structures at large lags 
by the French-hat filter or the broader peak without ripples by 
means of the Mexican hat. 

Because the field is periodic we can also apply the periodic 
continuation method here. The equivalent figure shows the same 
shape of the peak but more pronounced ripples with deeper min- 
ima on large scales for the French hat filter and a somewhat 
steeper decay for the Mexican hat. In this case the A-variance 
spectrum simply represents the Fourier transform of the filter 
function because the structure contains only a single Fourier 
component. For integer wavenumbers the periodic continuation 
is identical to the mirror continuation so that we obtain the same 
spectra except for a small discretisation error from a single pixel 
row which is treated different when mirroring or continuing peri- 
odically. In general the edge truncation always produces a strong 
and the mirror-continuation a weak smearing of the French-hat 
A-variance ripples at large lags relative to the periodic continua- 
tion. 

Corresponding results for the chess board map show almost 
the same curves as the sine wave field at all lags above the peak 
in the A-variance spectra. But they show a flatter rising part at 
smaller lags. The additional high-frequency contributions from 
the edges of the chess fields increase the A-variance at small 
lags leading to a slope shallower than 1 there. The situation is 
somewhat different for the circle map. The strong contributions 
from the "empty" area at large lags visible in Fig. |5] suppress 
differences from the different filter shapes in the spectrum. The 
spectra obtained from both filter types look very similar 

To quantify the influence of the selection of the filter shape 
and the diameter ratio on the detection of the dominant scale 
length we measure the sharpness of the A-variance peaks shown 
in Fig.|7] We use the width of the peak w, given as the logarithm 
of the ratio between upper and lower lag where the A-variance 
drops to 1/2 of the peak value, and its contrast relative to the 
values at lags which are either small or large compared to the 
peak position, i.e. relative to the first and the last point in the A- 
variance spectrum. The behaviour of both parameters is shown 
in Fig. [8] as a function of the filter diameter ratio v. The upper 
plot shows the logarithmic width of the peak; the lower plot the 
contrast of the peak relative to values at much larger and much 
smaller lags. From the figure it is obvious that we cannot achieve 
a minimum width and maximum contrasts simultaneously with 
the same filter, so that some balance has to be found. The con- 
trast with respect to smaller lags hardly varies with filter type 
and diameter ratio but the contrast relative to large lags is drasti- 
cally changed. The French-hat filter always produces a narrower 
peak but the decrease of the peak width towards lower v ratios 
is accompanied by a deterioration of the contrast with respect to 
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Fig. 8. Logarithmic FWHM of the A-variance peak and contrast 
between the peak value and the values at I = 0.014 and 1-1/3 
for the sine wave map analysed in Fig.|7] 



large lags. The Mexican-hat filter shows a continuous improve- 
ment of both quantities towards lower ratios but gives a some- 
what broader peak. The slight increase of the contrast with re- 
spect to the lower end of the spectrum towards lower v values 
for both filter types indicates that low diameter ratios result in a 
somewhat longer dynamic range below the dominant peak where 
the A-variance spectrum follows a power law. 

Comparing the different edge treatment approaches in cor- 
responding plots shows that the mirror-continuation always pro- 
duces about the same peak width and a somewhat better contrast 
than the filter truncation method because the latter introduces 
some flattening in the A-variance spectra at the largest lags. 

A very similar behaviour is also observed for the chess board 
structure. In contrast, the filled-circle map with its large-scale 
contributions shows a different behaviour demonstrated in Fig. 
|9] Only the contrast relative to the A-variance values at small 
lags shows the same small improvement for both filters towards 
low diameter ratios v as seen in Fig. [8] The other parameters 
behave qualitatively different. For the Mexican hat we find an 
optimum diameter ratio v ^ 1 .2 ... 1 .4 where the peak width has 
a minimum and the contrast with respect to large lags shows a 
maximum. For the French hat the contrast relative to large lags 
deteriorates over a much broader range at low v ratios and we 
observe an increased width of the peak there. Diameter ratios 
below V ^ 2.3 clearly reduce the sensitivity of the French-hat fil- 



Fig. 9. Logarithmic FWHM of the A-variance peak and contrast 
between the peak value and the values at I - 0.014 and Z = 1/3 
for the filled circle map discussed in Fig. |5] The widths mea- 
sured for the French hat at v < 1 .25 are lower limits because the 
spectrum did not fall below half of the peak value at the upper 
end of the spectrum. 



ter in this case. The increase of the peak width and the reduction 
of the contrast can be understood as resulting from the stronger 
pickup in the side lobes of the Fourier transformed French-hat 
filter which are closer to the main peak and stronger for lower v 
ratios. 

The corresponding figure computed using the mirror- 
continuation is very similar. The mirror-continuation always 
provides a somewhat better contrast. The eff'ect of contrast re- 
duction at large lags by the filter truncation is always higher 
for the French hat than for the Mexican hat and decreasing with 
growing distance of the filled circle from the map boundaries. 

The somewhat better contrast relative to small lags obtained 
in all cases with the French-hat filter indicates that this filter is al- 
ways more suitable to detect a power-law scaling behaviour over 
a wide dynamic range below a dominant scale. For the detection 
of the dominant scale the situation is less clear. In the sine wave 
and chess board maps the Mexican hat provides a better contrast 
but a larger peak width, for the circle map the Mexican hat at its 
optimum diameter ratio is only slightly worse than the French 
hat at its optimum ratio. Taking the general uncertainties from 
the French-hat ripples at large lags, however, it seems always 
preferable to use the Mexican-hat filter. 
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Comparing all results with respect to a clear indication of 
particular structure scales we find that the Mexican-hat filter 
with a diameter ratio v < 1.4 provides the best resolution. 
Diameter ratios v between 1 .4 and 1 .7 still produce a reasonably 
good sensitivity. The French-hat filter has its maximum sensitiv- 
ity for diameter ratios v between 2.3 and 2.5. Although it pro- 
duces A-variance peaks with a smaller width than the Mexican- 
hat filter, it produces in general ripple artifacts in the A-variance 
spectrum at large lags, lowering the overall contrast of the peak, 
so that it should be deferred relative to the Mexican-hat filter 
for the structure detection. For non-periodic structures covered 
by regular maps with rectangular boundaries, the edge treatment 
by mirror-continuation gives a somewhat better contrast than the 
filter truncation method but this difference is relatively small. 

4.3.2. Spectral index 

To judge the value of the different filters with respect to the re- 
trieval of the correct slope of an fBm structure or a sub-map from 
an fBm structure we consider two quantities: the dynamic range 
over which the A-variance slope can be reliably determined and 
the difference between the actually measured slope and the the- 
oretical value. 

The dynamic range of scales traceable with a filter of given 
shape and diameter ratio is constrained by the maximum filter 
size that can be used for the given map size. It is measured by 
the maximum effective filter length for which the contributions 
from those parts of the filter extending beyond a map boundary 
produce no noticeable distortion of the A-variance spectrum. In 
the case of the French-hat filter we find that the diameter of the 
annulus, i.e. / x v, may not exceed the size of the map for reliable 
A-variance values. For the Mexican hat, the parameter I x v must 
not exceed 2/3 of the map size. From the relations between the 
effective filter length and the filter size / obtained in Sect. 14.21 
we see that the dynamic range of effective lags available for a fit 
of the A-variance spectrum grows with decreasing diameter ratio 
V. Smaller v values increase the total range of lags where the A- 
variance can be determined without being dominated by edge 
effects. Moreover, Sect. l4.3.Tl shows that lower v ratios also tend 
to extend the dynamic range of scales below a structure peak 
where the A-variance follows a power law. Thus, low ratios also 
seem to be favourable for a reliable slope detection. 

To study the agreement of the measured A-variance slopes 
with the theoretically predicted index as a function of the filter 
shape we compute the spectra for fBm structures and sub-maps 
from fBm structures using the different filters. Fig.[TO]shows the 
spectra for a submap from an fBm with ^ = 3.0 computed with 
four different filters. First we notice the extension of the dynamic 
range for lower v ratios. None of the spectra gives an exact power 
law, but the French-hat filter with v = 3.0 and both Mexican-hat 
filters provide a reasonably good reproduction of the theoretical 
index a - 1.0. For low v ratios, the French hat tends to overes- 
timate the true spectral index. The Mexican hat results in some- 
what too low exponents. 

The corresponding spectra computed by the help of the 
mirror-continuation method result in a spectral index which is 
too low by about 0.1 compared to the theoretical value due to 
the flattening of the spectra at large lags as seen in Fig. [3] The 
Mexican hat and the French hat provide almost the same slopes. 

When applying the analysis to the sine wave field with k - \, 
i.e. the equivalent of an fBm with ^ = oo, the longer dynamic 
range traced by the French-hat filter results in a measured spec- 
tral index of the A-variance which is closer to the theoretical 
value of 4 than that obtained for the Mexican-hat filter. 
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Fig. 10. A-variance spectra computed for a sub-map from an 
fBm structure with ^ = 3. Filter truncation is used for the edge 
treatment, a denotes the power-law slope fitted over the full 
range. 



For a systematic investigation of the accuracy in the deter- 
mination of the A-variance slope as a function of the filter di- 
ameter ratio V we analyse sets of 30 different fBm structures 
and 30 submaps from fBms using the two basic filter shapes 
varying their diameter ratio v and the fBm spectral index ^. 
Fig. [TT] demonstrates the result for the French-hat filter ap- 
plied with the filter-truncation method. The range of spectral 
indices covers the typical indices in observations of inter stellar 
clouds jElmegreen & Scald 120041; iFalgarone etall |2004 . The 
strongest deviations of the measured exponents from the ex- 
pected values occur at low diameter ratios v. Hence, ratios below 
1.7 should not be used. fBms and submaps behave differently. 
For the fBms the strongest deviations from the theoretical spec- 
tral index occur at low spectral indices; for the fBm submaps 
they occur at high indices. This can be explained by selection 
effects. When selecting submaps from an fBm, there is a signif- 
icant scatter in the properties of the actually selected structure. 
This is visible as wider error bars. For high spectral indices we 
often find submaps which are dominated by some structure ex- 
tending beyond the submap boundaries. This tends to increase 
the measured slope. 

Other combinations of filter shape and edge treatment show 
somewhat different properties in details but the same general be- 
haviour as Fig. [TT] The mirror continuation method always un- 
derestimates the spectral index. With both filter types it produces 
at diameter ratios v > 2.0 slopes which are too low by 0.05-0.1. 
At lower ratios the difference grows up to 0.2. The Mexican-hat 
filter always shows a slightly stronger deviation from the theo- 
retical value than the French-hat filter. However, when systemat- 
ically correcting the slopes by a constant shift of Aor = 0.1 both 
filter types provide reliable spectral indices for v > 1 .4. With the 
filter truncation method both filters result in a good reproduction 
of the spectral index at diameter ratios v > 2.5. An acceptable 
sensitivity is still obtained for v > 1 .4 with the Mexican hat and 
for V > 2.0 with the French-hat filter. In this intermediate range, 
the measured average slope always deviates by less than 0. 1 from 
the theoretical value. 

One has to keep in mind, however, that most of the deviations 
discussed here fall below the size of the statistical error bars. For 
the fBm submaps, which are most representative for astronom- 
ical maps, they amount to Aa x +0.15 at a map size of 128^ 
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Fig. 11. Distribution of fitted A-variance exponents for fBm 
maps (upper plot) and submaps from fBm structures (lower plot) 
with ^ - 2.5, 3.0, and 3.5 as a function of the filter diameter ratio. 
The French-hat filter combined with the filter truncation method 
is used. Dotted lines indicate the expected A-variance exponents 
corresponding to the fBm spectral indices. The error bars indi- 
cate the scatter obtained for structures with the same parameters 
in terms of the standard deviation. 



and to Aa » +0.25 at a map size of 32^ (see Sect. 13. It . Except 
for very low diameter ratios, v < 1 .2, the use of the Mexican 
hat always provides a somewhat lower scattering of the mea- 
sured exponents than the French-hat filter. The edge treatment 
has practically no influence on the size of these scatter bars. 

Comparing the results from the scale detection and the repro- 
duction of the spectral index we find that each problem asks for 
a different optimum filter Whereas the scale detection favours 
the Mexican-hat filter with a diameter ratio v :$ 1 .4, the slope 
reproduction is equally well satisfied by both filter types for a 
diameter ratio v ^ 2.0. A reasonable compromise, providing a 
good sensitivity to either issue is the use of a Mexican-hat filter 
with a diameter ratio v » 1.5. 

The edge treatment by mirror continuation is somewhat 
favourable relative to the filter truncation in the slope detection, 
but requires that the average spectral index is corrected by a con- 
stant shift of Aa = 0.1. On the other hand the filter truncation 
method can be applied as well for maps with irregular bound- 
aries and shows a somewhat lower statistical scattering within 



the studied samples. Thus, the filter truncation provides the most 
reliable parameters in a single run for any data set without the 
need for additional corrections. 



5. Conclusions 

The A-variance analysis was previously established as a general 
tool to study the scaling behaviour of interstellar cloud struc- 
ture. The main advantage of the A-variance method compared 
to the computation of the power spectrum is its robustness with 
respect to angular variations, singular distortions, gridding and 
finite map size effects. 

We propose two essential improvements of the A-variance 
analysis. The first one is the use of a weighting function for each 
pixel in the map. This weighting function allows us to study data 
sets with a variable data reliability across the map and to simul- 
taneously solve boundary problems even for maps with irregular 
boundaries. 

Maps with a variable data reliability are eventually obtained 
in most observations, either due to a local or a temporal vari- 
ability of the detector sensitivity or the atmosphere or due to 
different integration times spent for different points of a map. 
The ordinary A-variance analysis as well as the power spectrum 
fail to take the resulting effects into account. By applying the im- 
proved A-variance analysis to noisy data we find that only the use 
of a significance function to weight the different data points al- 
lows us to distinguish the influence of variable noise from actual 
small-scale structure in the maps. In the case of statistical un- 
certainties or sensitivity changes the weighting function is best 
provided by the inverse RMS in the data points. The weight- 
ing function allows us to considerably extend the dynamic range 
within which the intrinsic scaling behaviour of an observed as- 
tronomical structure can be measured (Appx.lAli. 

In the treatment of map boundaries the use of a weighting 
function allows us to use computational methods like the fast 
Fourier transform to compute the A-variance spectrum by ex- 
tending a measured map by points with zero significance. This 
is mathem atically equi valent to the truncation of the filter as pro- 
posed bv iBensch et al J (1200 1 ) but has the advantages of the fast 
computation and the definition of a smooth filter shape in Fourier 
space not influenced by gridding effects in ordinary space. The 
virtual filter truncation is the only approach to analyse maps 
which are only sparsely filled by significant values. For rectan- 
gular maps a periodic continuation by mirroring can be used as 
well to solve the boundary problem. Mirror-continuation will al- 
ways underestimate the spectral index by about 0. 1 . This could 
be easily taken into account. However, the mirror-continuation 
is less flexible than the filter truncation which works on irregular 
maps as well although resulting in slightly wider statistical error 
bars. 

The second improvement of the A-variance analysis is its op- 
timisation with respect to the shape of the wavelet used to filter 
the observed maps. We have computed the effective filter length 
as a function of the shape and compared the peak positions of 
resulting A-variance spectra with characteristic structure sizes in 
the test data sets. We find that the peak positions always falls 
10-20 % below the maximum structure size. Taking this system- 
atic offset into account we can calibrate the spatial resolution 
of the A-variance analysis to about +5 %. Unfortunately, it is 
not possible to define a single optimum wavelet for all purposes 
because different wavelets show different qualities in the detec- 
tion of the characteristic structure scaling behaviour The best 
choice for an exact measurement of the power spectral slope 
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Fig.A.l. A-variance spectra computed from an fBm structure 
with a spectral index ^ = 3.1 and a superimposed circular pattern 
of white noise. The solid line shows the A-variance spectrum of 
the original fBm, the dotted line the spectrum obtained from the 
direct analysis of the noisy map, and the dashed line the spec- 
trum computed with a reliability weighting given by the inverse 
noise RMS. 



are wavelets with a high ratio between the diameter of the an- 
nulus and the core of the filter. Here, the French-hat and the 
Mexican-hat filter are equally well suited. For the detection of 
pronounced size scales the Mexican-hat filter and low diameter 
ratios are preferred. A good compromise between the different 
requirements is the Mexican-hat filter with a diameter ratio of 
1 .5 always providing a A-variance spectrum with approximately 
the correct slope and without missing any special spectral fea- 
ture. 

We provide an easy-to-use IDL widget program implement- 
ing the two-dimensional A-variance analysis as described here 
implementing the different filters and edge-treatment methods 
for the analysis of arbitrary maps in FITS formafl 
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Appendix A: The impact of a varying data reliability 

As a combined test of all extensions we use the results on the op- 
timum filter shape with the weighting function and return to the 
analysis of data with a varying reliability within the map. To test 
how the improved A-variance analysis recovers the properties of 
an original structure from measurements influenced by a vary- 
ing noise level, we create maps where white noise with a spatial 
pattern of different noise amplitudes was added. We use combi- 
nations of the different spatial structures discussed in Sect. 12.21 
for the structure to be measured and for the spatial distribution 
of a noise level superimposed to the data. 

Fig. |A. II shows one example of a resulting A-variance spec- 
trum for an fBm structure with a spectral index ^ = 3.1 where 
a noise pattern given by the filled circle described in Sect. 12.2.21 
and d - 213 was added. The average signal-to-noise ratio, de- 
fined as the ratio between the maximum in the fBm structure 
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and the noise RMS, is 1 and the variation between the noise lev- 
els inside and outside of the circle is a factor 9. This example 
may represent the situation of an observed map where the in- 
ner part is covered by many integrations, so that it shows a high 
signal-to-noise ratio, whereas the outer part is observed with few 
integrations leading to a higher noise level. 

The solid line represents the A-variance spectrum of the orig- 
inal fBm structure. The dotted line is the spectrum that is ob- 
tained by the direct analysis of the noisy map without any re- 
liability weighting. Because there is no correlation in the noise 
between neighbouring data points, the added noise contributes 
to the A-variance spectrum only on small scales with a decay 
proportional to l"^ towards larger lags. Due to the relatively high 
average noise level, the A-variance spectrum is dominated by 
the noise contribution up to lags of about 0.2. The original fBm 
spectrum is only matched within a very narrow scale range at the 
largest lags. 

The dashed line shows the improvement that is obtained by 
using the knowledge on the noise level in terms of a weighting 
function w{r) inversely proportional to the local noise RMS. Due 
to the relative suppression of contributions from the outer noisy 
parts of the map the original fBm spectrum is recovered over a 
much broader range of scales. We find, however, a small distor- 
tion at the data point for the largest lag. This can be interpreted as 
the effect of a slight "cross talk" from the weighting function to 
the measured structure. The obvious strong improvement of the 
recovery of the original A-variance spectrum from the noisy data 
is thus achieved at the cost of a slightly reduced data reliability 
on the characteristic scales of the weighting function. 

To study this effect more systematically we perform a num- 
ber of parameter studies combining the different structures with 
varying noise patterns, varying noise levels and varying noise 
dynamic ranges. In the resulting A-variance spectra we com- 
puted the scale range over which the spectrum agrees with the 
spectrum of the original structure within 10%. The length of 
this range, which allows a reliable derivation of the true scaling 
behaviour, is a measure for the quality of the structure recovery. 
To test the influence of selection effects we repeat each computa- 
tion for a number of different initialisers for fBm structures and 
for the noise fields, so that we arrive at 30-80 computations for 
each parameter set providing a statistically significant sample. 

The result of such a parameter scan is demonstrated in Fig. 
IA.2I In this example, the original structure is an fBm structure 
with f = 3.0 and the superimposed noise amplitude follows a 
chess board structure (see Sect. l2.2.TI ) with four fields, i.e. a field 
length of half the map size. By changing the noise amplitude 
ratio between the different fields while preserving the average 
noise amplitude this plot is suited to study the actual effect of the 
noise variation and the corresponding correction by a weighting 
function in the A-variance analysis. 

In the analysis without weighting function, the dynamic 
scale range within which the true spectrum can be fitted always 
covers about a factor 9, independent of the variation in the noise 
level between the four fields. The noise correction by the weight- 
ing function can extend this range up to an average factor of 25 
in the case of high variation levels. The error bars, indicating the 
minimum and maximum ranges detected in the sample, show, 
however, that there is a considerable spread in the range length 
over which the fit is reliable. The matching range is always in- 
creased compared to the A-variance analysis without weighting 
function, but the actual magnitude of this increase can consid- 
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Fig.A.2. Scale range (ratio between maximum and minimum 
lag) within which the derived noisy spectrum agrees within 10 % 
with the original spectrum when the A-variance spectrum is 
computed either with or without a weighting function given by 
the inverse noise RMS. The original structure is given by an 
fBm structure with ^ = 3.0 and the noise pattern by a chess 
board structure with four fields and a signal-to-noise ratio of 
5.0. On the abscissa we have changed the amplitude of the noise 
level variation between the different fields of the chess board. 
The error bars show the absolute minimum and maximum of the 
matching range size found in a sample of 30 different fBms. 



coverage of the main features of a structure by low-noise ob- 
servations. Future studies should show whether the information 
from the A-variance spectrum of the map of weights can be used 
to further improve the outcome. 
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erably vary0. The example represents, however, a kind of worst 
case scenario, because the noise variation in this pattern virtu- 
ally cuts out two large pieces from the map which may contain 
main elements of the original structure. This is not expected for 
real observations where the astronomer would hardly select a 
field avoiding the main object of interest. For a chess board noise 
with half the cell size, the error bars for the distribution of fitting 
ranges are already reduced by almost a factor two. The example 
is nevertheless instructive because it shows all the effects that 
we encounter in the parameter study with different intensity and 
noise structures. We always find the increase of the fitting range 
but also the wider spread of the ranges within the sample studied. 
When using fBms for the noise distribution with spectral indices 
( around 4 or higher, we find as well a wide spread of the factors 
by which the length of the fitting range is increased in the new 
A-variance analysis. For all lower indices the error bars shrink 
in the same way as for the chess board with smaller cell size. 
In general, the extension of the matching range by the use of a 
weighting function is most significant for strongly fragmented 
noise maps and noise maps well adapted to the dominant parts 
of the actual structure. The latter case is usually given in astro- 
nomical observations. 

We conclude that the introduction of a weighting function 
given by the inverse noise of the data into the A-variance analysis 
always extends the spatial range over which the original scaling 
behaviour can be recovered. The amount of this improvement 
depends on the strength of the noise variation across the map 
and the coverage of the observed structure by regions of low 
noise. The range within which the measured A-variance spec- 
trum agrees with the original spectrum can be extended by more 
than a factor three for high levels of noise variations and a good 



^ The upper limit at about 30 corresponds to the whole spectrum. 
Thus a further extension is not possible here. 



